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We consider domain walls between regions of superfluid ^He-B in which one 
■ component of the order parameter has the opposite sign in the two regions 

^ ■ far from one another. We report calculations of the order parameter profile 

O I and the free energy for two types of domain wall, and discuss how these 

structures are relevant to superfluid ^He confined between two surfaces. 
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\ The order parameter of superfluid He is suppressed near an impurity 

' or in the vicinity of walls and interfaces .^^^'^ This happens because the 

order parameter breaks the rotational symmetry of the normal phase, and 
', because particle-hole coherence is destroyed when quasiparticles, which move 

in the presence of the order parameter, are scattered from one point on the 
Fermi surface to another. Here we consider a related problem in which the 
order parameter varies strongly in space, so that quasiparticles move in this 
\ spatially varying order parameter even when they propagate along straight 

trajectories. Consider an order parameter that varies along some axis and 
Q . has one bulk solution on the far left, and another possible bulk solution 

on far right. The region separating these two degenerate bulk solutions is 
the domain wall. Since the order parameter is different on the two sides 
the superfluid will be strongly deformed in the vicinity of the domain wall. 
5^ \ Although domain walls are not energetically favorable in bulk ^He, we will 

argue that inhomogeneous structures formed from domain walls may be the 
favorable solutions for ^He in confined geometries. 

To determine the order parameter and to calculate the corresponding 
free energy of a domain wall we use the quasiclassical theory.!^ We solve the 
Eilenberger transport equation for the quasiclassical propagator, ^(p, R; 
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in the Matsubara formalism, 

[iEmTs -A,g]+ivf'Vg = 0, (1) 

where £m = (2m+l)7rr, A is the 4x4 matrix order parameter and vj = Vfp 
is the Fermi velocity corresponding to the point pj on the Fermi surface. 
In addition, the physical solutions of the quasiclassical propagator satisfy 
Eilenberger's normalization condition, g'^ = — vr^l. It is convenient to express 
the propagator and order parameter in particle- hole and spin space as, 

g f • (iaay) \ ^ _ f {i(ray) ■ A 



(2) 

The vector in spin space, A(p,R), depends on the spatial position, R, and 
the direction of the momentum, p. For pure p-wave pairing the order pa- 
rameter is traditionally written in terms of the spin (index a) and orbital 
(index i) matrix, Aai(R), 

3 

Aa{p,K)=Y,Aa^{^)pi. (3) 

i=l 

The gap equations for the order parameter components, 
^«.(R)ln^ = 3rg [^p.x 



e,«>0 



X ( f (p, R; + f'(p, R; - 27rM^ ) , (4) 



must be solved self-consistently with the transport equation. 

Let us choose the x-axis to be the axis along which one of the compo- 
nents of the order parameter is changing. The domain wall is in the yz-plane 
centered at x = 0. We assume the order parameter to be translationally in- 
variant in the yz-plane. Reflection symmetries under y —y and z — > —z 
allow us to write the order parameter in a simple diagonal form, 
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= ylj^y I = I Ai II I , (5) 
\ 

or equivalently in vector form. 



(6) 
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where the symbol || (_L) refers to momentum directions parallel (perpendic- 
ular) to the domain wall. 

We consider two basic configurations of domain wall. The bulk phase 
of superfluid '^He at zero pressure is the B-phase with the order parame- 
ter Ao(p) = (AqPj^, AoPy, Other possible bulk solutions differ from 
this one by a sign change of one of the components. We consider the first 
domain wall configuration to be one with an order parameter A.{p)l = 
(-AoPx, AoPy, Aop^) on the far left, and Mp)r = {+Aop^, AoPy, Aqp^) 
on the far right. We refer to this configuration as the "perpendicular" do- 
main wall because it is A^ that changes sign across the domain wall. This 
configuration is equivalent to the problem of determining the order param- 
eter structure near a specularly reflecting surface positioned at x = 0. In 
the case of a specular surface a quasiparticle with momentum (pxiPyiVz) is 
scattered to a state with momentum {—Px,Py,Pz)- As a result these spec- 
ularly reflected quasiparticles move in the same order parameter field as 
quasiparticles moving along straight trajectories through a perpendicular do- 
main wall. The second domain wall configuration we consider has A(p)i = 
(AoPx, AoPy, -Aop^) on the left and A(p)r = (AqPx, AqPj/, +Aop^) on the 
right side, and we refer to this structure as a "parallel" domain wall. This 
structure is not related any order parameter structure defined by a refiecting 
surface. 

The self-consistent solutions for the two order parameter configurations 
are shown in Fig. ^ We see when Ay changes sign (right panel) the domain 
wall is approximately twice as narrow as than the domain wall in which 
A_L changes sign. This means that the deformation of the superfiuid order 
parameter is weaker for the case of the "parallel" domain wall. The differ- 
ence in the order parameter structures between these two domain walls is 
qualitatively explained by difference in the number of trajectories that are 
pairbreaking. In the "perpendicular" configuration the pairbreaking trajec- 
tories are predominantly perpendicular to the domain wall. In the case of the 
"parallel" domain wall trajectories that are nearly parallel to domain wall 
are most pairbreaking. However, these trajectories only weakly connect the 
two sides of the domain wall. As a result the profile of the order parameter 
is narrower than that of the perpendicular domain wall. 

We explore this feature further by computing the free energy for each 
given order parameter profile. The free energy is calculated by introducing 
the auxiliary propagator 'gx that satisfies the transport equation. 



i£mT3 - XA,gx 



+ iwf ■ Vg\ = , 



(7) 



with < A < 1. The free energy in the weak coupling limit can be expressed 
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Fig. 1. The order parameter profile across the domain wall (A is in units 
of 27rfcsTc, ^0 = hvf /2TTkBTc). Left panel: change in A_l across the wall, 
Axx{—x) = —Axx{x). In this case A;^|| = A211 = Ay. Right panel: A211 
changes sign, Azzi—x) = —Azz{x), and A^y is slightly different from A_l. 
For comparison the thin dotted line is A^ from the left panel. The parallel 
domain wall is thinner than the "perpendicular" domain wall by a factor of 
about 1.8. 




as an integral over the variable coupling constant, Xp^ 

^^ = \f^ Sp' I A - ^5) } ^ I d'R AF(R) . (8) 

We define the free energy density, AF(R), from the definition of the trace, 
Sp', which stands for summation and integration over all indices and vari- 
ables, e.g. 

Sp'{a?} ^ I d^ni^TY^Nf I ^Tr4(A(p,R)5(p,R;e„.))| • (9) 

The free energy densities for the two domain walls are shown in Fig. 
121 The change in A211 takes place over a shorter distance than that of A_|_. 
Thus, the former configuration has larger gradient energy. Nevertheless, this 
is compensated by the spatial extent of "perpendicular" domain wall with 
the result that the net energy cost of the "parallel" domain wall is less than 
that of the "perpendicular" domain wall. This result has implications for the 
structure of the order parameter and the stability of new phases of superfluid 
^He in confined geometry. 
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Fig. 2. Reduction of the condensation energy near a domain wall [in units 
Nf{2TTkBTc)'^]. The magnitude of the free energy density for the "parallel" 
configuration is suppressed more around x = due to the larger gradient 
energy terms, but the integrated suppression is less than that for the "per- 
pendicular" configuration. 

We demonstrate this point by the following argument. Fig. El shows a 
slab of ^He-B and two types of quasiparticle trajectories. Trajectory 1 under- 
goes a reflection from a boundary, while trajectory 2 is a straight trajectory 
in the region of interest. In case (a) shown in Fig. |31 the order parame- 
ter component A^z has the same sign to the left and right of the vertical 
line. Thus, quasiparticles moving on trajectory 1 propagate in a changing 
order parameter field, —^z, which is equivalent to the "perpendicu- 

lar" domain wall. However, quasiparticles moving along trajectory 2 do not 
experience a change in the sign of the order parameter. 

Case (b) shown in Fig. |S1 is quite different. In this case we assume a 
domain wall in which Azz changes sign on the two sides. Now quasiparticles 
moving on trajectory 1 do not experience a sign change because the change 
in the momentum, pz —Pz, is compensated by the sign change in Azz- 
Thus, there is no contribution to the order parameter deformation from this 
trajectory. However, this reduction in surface pairbreaking occurs at the cost 
of pairbreaking introduced trajectories of the second type. Quasiparticles 
moving on trajectory 2 propagate through an order parameter field with a 
sign change. But since these type of trajectories lead to a domain wall of the 
"parallel" type, the domain wall structure may be favorable to the uniform 
solution in the thin film shown in case (a). 

This observation about the possible stability of domain wall structures 
depends on the thickness of the slab or film since for very thick films it is 
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Fig. 3. Pairbreaking trajectories with and without a domain wall in the 
■^He-B film, a) no domain wall, b) Azz{x) changes sign from left to right. 

clearly unfavorable to have domain walls in the plane of the film. How- 
ever, for thinner films or slabs inhomogeneous structures related to parallel 
domains are possible. A full quasiclassical calculation of the structure and 
free energy for thin films of superfluid ^He confirms this observation based 
on the domain wall calculations. We have found that associated with the 
domain wall structures new order parameter components develop, and that 
there are periodic structures formed from domain walls along the plane of 
the film that are energetically favorable compared to an order parameter 
that is translationally invariant in the plane of the film. We find a range of 
film thickness, lO^o ~ I^^Qj for which these inhomogeneous states are favored. 
The details of this calculation will be published elsewhere. 

We thank Tomas Lofwander and Erhai Zhao for useful comments on 
this work. 
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